Charge Frustration Effects in Capacitively Coupled Two-Dimensional 

Josephson-Junction Arrays 
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We investigate the quantum phase transitions in two capacitively coupled two-dimensional 
Josephson-junction arrays with charge frustration. The system is mapped onto the S — 1 and 
S — 1/2 anisotropic Heisenberg antiferromagnets near the particle-hole symmetry line and near 
the maximal-frustration line, respectively, which are in turn argued to be effectively described by a 
single quantum phase model. Based on the resulting model, it is suggested that near the maximal 
frustration line the system may undergo a quantum phase transition from the charge-density wave 
to the super-solid phase, which displays both diagonal and off-diagonal long-range order. 



PACS numbers: 74.50. +r, 67.40. Db, 73.23.Hk 

In recent years, various types of cotunneling transport 
have been of great interest in ultrasmall tunnel junct. 
tions, which exhibit strong Coulomb blockade effect .El 
In particular, cotunneling of the electron-hole pairs in 
two capacitively coupled one-dimensional (ID) arrays 
of small metallic junctions has been proposed theoreti- 
callyo and demonstrated experimentally!] revealing the 
remarkable effects of the cucsent mirror. More recently, 
in capacitively coupled IDqB or two-dimensionaB (2D) 
Josephson-junction arrays (JJAs), the cotunneling of 
particle-hole pairs (with the particle and hole standing 
for the excess and deficit Cooper pair, respectively) has 
been proposed even to drive the quantum phase transi- 
tion from superconductor to insulator (SI) at zero tem- 
perature. Here the particle-hole symmetry of the sys- 
tem may be broken by, e.g., the gate voltage applied be- 
tween the array and the substrate. The resulting charge 
frustration is expected to affect the phase transition of 
the system in a crucial way. For example, when the 
particle-hole symmetry is broken maximally, the trans- 
port is governed by the cotunneling of the particle-void 
pairs (with the void denoting the absence of an excess 
or deficit Cooper pair) and the different nature of the 
associated .phase transition has been pointed out in one 
dimension. □ On the other hand, existing studies of cou- 
pled 2D arrays with charge frustration have concentrated 
upon the charge- vortex duality,0 without appreciable at- 
tention to the phase transitions. n 

In this paper, we extend the previous workQ on two 
capacitively coupled 2D arrays of ultrasmall Josephson 
junctions to investigate the charge-frustration effects on 
the quantum phase transitions. In a manner similar 
to that of Ref. |J, we map the system to the S = 1 
anisotropic Heisenberg antiferromagnet near the particle- 
hole symmetry lines and to the S — 1/2 one near the 
maximal- frustration lines. It is then argued that the two 
spin models can in effect be incorporated into a single 2D 
quantum phase model with the effective self-capacitance 
given by the coupling capacitance of the original two- 
array system and the junction capacitance by the intra- 
array junction capacitance. The resulting model indi- 
cates that near the maximal frustration line the system 
may exhibit a quantum phase transition from the charge- 



density wave (CDW) to the super-solid (SS) phase. In 
the SS state, the system possesses both diagonal and off- 
diagonal long-range order (DLRO and ODLRO): Namely, 
both the density-correlation of charges and the phase- 
correlation of superconducting order parameters remain 
finite as the distance grows arbitrarily large.tm 

The system of coupled 2D square arrays, shown 
schematically in Fig. [j], is described by the Hamiltonian 

H = 2e 2 J2 K(r)-%]C„;(r,r')K(rVR 9 ] 

e,£';r,r' 

- Ej ^2 cos[<fo(r) - <^(r+e M )] 



H c + Hj, 
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where the number ne(r) of the Cooper pairs and the 
phase 4>i(r) of the superconducting order parameter 
at site r on the £th array (I = 1,2) are quantum- 
mechanically conjugate variables: [ng(r), <^(r')] — 
iSwS rr '. The Josephson coupling between neighboring 
sites r and r+e M (with e M being the unit vector in the 
direction /i = x, y) in each array is characterized by 
the coupling energy Ej, whereas the external charge 
n g = CoV g /2e induced on each island by the applied gate 
voltage V g breaks the particle-hole symmetry of the sys- 
tem, introducing charge frustration. The two arrays are 
coupled through the capacitance Cj between two grains 
at the same position r on the two arrays. (Note the diffcrp. 
ence from the Josephson coupled multi-layered system,EI 
where Cooper-pair tunneling between layers is allowed.) 
The capacitance matrix C characterizingthe charging en- 
ergy part He of the Hamiltonian in Eq. (yj) can be written 
in the block form: 



,(t,t>) = C(t,t> 
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1 



>Cr 



1 -1 
-1 1 



(2) 



where C(r, r') is the usual intra-array capacitance matrix 
C{t,t') = C 8 ttI 



Ci [ 2S * 



Jr,r'+e M J 



(3) 
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with Co and C\ being the self- and junction capacitance, 
respectively. Although it is not essential in the subse- 
quent discussion as long as the interaction range is fi- 
nite, we assume for simplicity that Ci/C'o < 1, keeping 
only the on-site and the nearest-neighbor interactions be- 
tween the charges. We also define charging energy scales 
E a = e 2 /2C* , E x = e 2 /2Ci, and Ei = e 2 /2C/, associ- 
ated with the corresponding capacitances. 

In the regime of concern in this paper, Cr 3> Co (> 
Ci), i.e., Ej <C E (< Ei), the charging energy part of 
the Hamiltonian in Eq. (pi) can be written conveniently 
as the sum He = H^+H^ with each component defined 
to be 

tf+^{/o$> + (r)-2n g ] 2 

r 

+ UiJ2[n+( r ) - 2n s ][n+(r+e M ) - 2n g ] (4) 
H c = V ]>>-(r)] 2 + H^n-Wiufr+e,), 

r r,j.L 

where we have defined new charge variables n±(r) = 
rii(r) ± 712(1") and the interaction strengths are given 
by U ~ 2E , V Q ~ E T , U x ~ 4(C 1 /C )B 0) and Vi ~ 
{C\/Ci)Ej. This representation of the charging energy 
part allows us to distinguish clearly the two interest- 
ing regions from each other: near the particle-hole sym- 
metry line n g — and near the maximal-frustration line 
n g = 1/2, as one can see from the energy spectra of 
He displayed in Figs. |2| and || for the two regimes, re- 
spectively (recall that Uq ^> Vq). (Since the system is 
invariant with n g — > n g + 1 , we need to consider only the 
range < n g < 1). As pointed out for two coupled ID 
arrays in Ref. |I| there follow the remarkable properties of 
the spectrum of He in each regime: Near the maximal- 
frustration line, the charge configurations which do not 
satisfy the condition n+(r) = 1 (for all r) have a huge 
excitation gap of the order of Eq. (Note that we are 
interested in the parameter regime £7, Ej « J?o.) Fur- 
thermore, the ground states of He, separated from the 
excited states by the gap of the order of Ej, have two- 
fold degeneracy for each r, corresponding to n_ (r) = ±1. 
This degeneracy is lifted by the Josephson coupling term 
Hj in Eq. (|l|) as Ej is turned on. As a result, it is con- 
venient in this case to work within the reduced Hilbert 
space £d, where 7i+(r) = 1 and n_(r) = ±1 for each 
r. Near the particle-hole symmetry line, on the other 
hand, the low-energy charge configuration should satisfy 
the condition 71+ (r) = for all r. Unlike the former case, 
the ground state of He is non-degenerate and forms a 
Mott insulator characterized by ni(r) = n 2 (r) = for 
all r. As Ej is turned on, the ground state of He is 
mixed with the states with n_(r) = ±2. Accordingly, 
the relevant reduced Hilbert space is given by £ s , where 
ra+ (r) = and 7i_(r) = 0,±2 for all r. 

Accordingly, it is instructive to project the Hamilto- 
nian in Eq. (|J) onto £ s (£ d ) for n g < 1/4 (for \n g -l/2\ < 



1/4); this results in the effective Hamiltonian, up to the 
second order in Ej/E , 



H, 



'41 



P 



H + H, 



1 - P 



E — H, 



-H, 



c 



P. 



(5) 



where P is the projection operator .0 Explicit implemen- 
tation of the projection near the particle-hole symmetry 
line can be achieved by first noting the correspondence 
between the charge picture of the original model and the 
pseudo-spin (S = 1) picture in the reduced Hilbert space 



S'(t) 



p ni(r)-n 2 (r) p 



S+(r) = v^Pe-^^l - F)e +,fc(r) P 
S~(r) = V2Pe- #2(r) (l - P)e+* (r) P. 



(6) 



In particular, the spin-flip operators S + and S~ manifest 
the second-order cotunneling process of the particle-hole 
pairs via an intermediate virtual state, as depicted in 
Fig. |](a), and mix the energy levels with unpaired par- 
ticles or holes of energy Uq and those with particle-hole 
pairs of energy 4Vo (see Fig. 0). It then follows that the 
effective Hamiltonian in Eq. ^) takes the form 

r 

- ^E^+M^r+e,) + S-(r)5+(r+e^)} , (7) 

which describes the spin-1 2D XXZ antiferromagnet.Eil 
Here the exchange interaction and the anisotropy ratio 
are given by 



7 E J i 1 
J ~ — — and 71 ~ — r, 

AE K 2 ' 



(8) 



respectively, with the dimensionless coupling constant 
K = y/ E 2 , 1 32EjEq . Near the maximal-frustration line, 
on the other hand, the effective Hamiltonian reduces to 
that for a spin- 1/2 2D XXZ antiferromagnet 



H 



^ /2 =7iJE^(r)^(r+e, 



J J EE {S + (r)S-(r+e,) + S~ (r)S+ (r+e„)} , (9) 
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with the exchange interaction and the anisotropic ratio 
given by 



3 ~ 2Gj K 2 ' 



(10) 



respectively. In this case, the definitions of the pseudo- 
spin operators in terms of the phase and charge operators 
are also different slightly from those in Eq. (o): 
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S z (r) - P ni ( r )~ n2 ( r ) P 

S+{r) = Pe- i,h{r \l - p) e +^ r )p 
S~(r) = p e -*( r )(l - p) e +*( r )p. 



(11) 



Such spin-flip operators are associated with the cotun- 
neling of the particle- void pairs as displayed in Fig. [|(b). 

In two dimensions, unlike the ID case, neither of 
the two (spin-1 and spin-1/2) XXZ antiferromagnets de- 
scribed by Eqs. (|7]) and (|^) allow exact solutions. The 
simple mean-field theory based on the Ginzburg-Landau 
approach^ indicates that the spin-1 XXZ antiferromag- 
net may exhibit a zero-temperature phase transition from 
the XY-like phase to the spin-1 Ising-like phase at K ~ 1 
or 71 ~ 1. In the charge picture, the XY-like phase cor- 
responds to the superconducting (SC) state displaying 
ODLRO while the spin-1 Ising-like phase characterized 
by S z (r) = describes the Mott insulator (MI) state with. 
DLRO. On the other handj-^nean-field-like approaches£3 
and numerical approachescJ to the spin-1/2 XXZ an- 
tiferromagnet suggest a zero-temperature phase transi- 
tion from the XY-like phase to the spin-1/2 Ising-like 
phase with staggered magnetization at K ~ yj Ci/2Cj 
or 71 ~ 1, corresponding to the SC state and the CDW 
state, respectively. 

Therefore, for the present, the projection of the Hamil- 
tonian to get the effective spin model does not provide 
us with direct information about the phase transitions. 
Remarkably, however, the spin models, given by Eqs. (]?]) 
and (^|) in the two regimes, can be obtained from a single 
2D quantum phase model (QPM), via appropriate pro- 
jections. This strongly indicates that both regimes can 
be described by the Hamiltonian for the QPM: 

H QPM = 2e 2 ^[n(r) - n„]C- 1 (r,r / )[n(r / ) - n g ] 



3_ 

4Pn 



5^oos[0(r)-0(r+e„)], (12) 



where the effective Josephson-coupling energy Ej/AEq 
is much reduced compared with the original intra-array 
value Ej, and the effective capacitance matrix reads 



C(r,r') = dS r 



(13) 



Note that the self-capacitance is given by Cj instead of 
the original value Co while the junction capacitance is 
given by Ci/2. The value of charge frustration fi g is re- 
lated with that (rig) of the original model given by Eq. (|l|) 
in the following way: At the symmetry line (n g = 0) and 
the maximal-frustration line (n g = 1/2) of the original 



model, we have n„ 



Near those lines, however, the 



value of fig is rather insensitive to that of n g : Namely, 
fig remains close to zero and to 1/2 in the rather large 



ranges around n g = and 1/2, respectively, changing its 
value sharply near n g w 1/4. Accordingly, the QPM in 
Eq. ( |l2] ) is either near the symmetry line (fi g ~ 0) or 
near the maximal-frustration line (n g as 1/2) except for 
the more or less narrow range around n g — 1/4. 

The reduction of the above QPM to the spin-1 and 
the spin-1/2 XXZ models via appropriate projections 
can be recognized as follows: In the case fi g ss 0, the 
charging energy reaches its minimum at n(r) = 0. This 
ground state becomes mixed with the states n(r) = ±1, 
as the Josephson coupling is turned on. On the other 
hand, for fi g rs 1/2, the minimum of the charging energy 
arises at n(r) = 0, 1, yielding two- fold degenerate ground 
states. These situations are essentially the same as those 
of the original model with charge frustration n g . We thus 
project the QPM onto the spaces £ s = {n(r) = 0, ±1} 
and £d = {n(r) = 0, 1} with the psuedo-spin operators 
redefined as 



S z {r) ee Pn{r)P 
S+{r) ee VaPe-^Wp 
S~(r) ee V2Pe^ (r) P 



(14) 



in space £ s and 



S z (r) = Pn(r)P- 1/2 
S+{r) ee Pe-^Wp 
S-(r) ee Pe l0(r) P 



(15) 



in space £d\ these projections reproduce, in the zeroth 
order of Ej/Ei, both the spin-1 and the spin-1/2 XXZ 
models in Eqs. (0) and @ for n g < 1/4 and \n g -l/2\ < 
1/4, respectively. 

As we proceed to higher orders, the projection of the 
single-layer QPM in Eq. ( [l2"| ) in general yields the coeffi- 
cients of the nth-order terms (Eq/ Eip n ~^ times larger 
than those in the projection of the original model in 
Eq. (|l]). In spite of such discrepancy in numerical coef- 
ficients, the two projections (of the original model and 
of the single-layer QPM) onto their own spin models 
should bring about quite similar structures. For exam- 
ple, mixing of the energy levels in £ s with those satisfying 
n_(r) = ±4, ±6, • • • (but still keeping n+(r) = 0) always 
occurs via the virtual states with energies of the order of 
Eq. Consequently, at least in the two regimes of concern 
here, it is not irrelevant to consider the single-layer QPM 
in Eq. ([IJ) as an effective model for the original system. 
Quite naturally, the deviation of the QPM in Eq. ( |T^ ) 
from the original model increases with Ej. 

The 2D QPM has been studied extensively in recent 
years (see, e.g., Ref. ||). Remarkably, for \n g — 1/2] <C 1/4, 
it was suggested that there may exist an unusual SS phase 
with both the DLRO and ODLRO, i.e., the coexistence 
of the crystalline charge ordering together with supercon- 
ductivity. The existence of the SS phase conflicts with 
the prediction of a direct transition from the CDW to 
the SC based on the spin-1/2 XXZ model in Eq. (§), 
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but such conflict also appears when one simply truru 
cates the effects of higher energy levels in the QPM.Ej 
These arguments finally yield the schematic phase dia- 
gram shown in Fig. where the thick solid lines repre- 
sent the phase boundaries of the SI transitions, separat- 
ing the SC from the MI (near the symmetry line) or from 
the CDW (near the maximal-frustration line depicted 
by the dashed-dotted line). Note that these boundaries 
(near the two lines) change rather gradually as n g is var- 
ied, which reflects that near the two lines the effective 
charge frustration h g in the QPM is insensitive to the 
original charge frustration n g . The dashed lines in Fig. ^ 
represent the somewhat speculative boundaries discussed 
above; here the region occupied by the SS phase might be 
small because in the QPM the self-capacitance is .much 
larger than the junction capacitance (Cj 3> Ci/2).0 

In conclusion, we have investigated the quantum phase 
transitions in two capacitively coupled two-dimensional 
Josephson-junction arrays with charge frustration. The 
system has been mapped into the S = 1 and the S = 1/2 
anisotropic XXZ antiferromagnets near the particle-hole 
symmetry line and the maximal-frustration line, respec- 
tively. We have then argued that the two spin models in 
effect can be incorporated into a single quantum phase 
model. Based on the resulting model, it has been sug- 
gested that near the maximal frustration line the system 
may exhibit a quantum phase transition from the charge- 
density wave to the super-solid phase, displaying both 
diagonal and off-diagonal long-range order. 
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FIG. 1. Schematic diagram of the two coupled 2D arrays. 
Each of the upper and lower arrays represents a 2D array 
composing the system. 
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FIG. 3. Energy levels and corresponding charge configura- 
tions near the maximal-frustration line. It should be noticed 
that the ground state is two-fold degenerate per site. 




(a) 





FIG. 2. Energy levels of He in Eqs. (^) and corresponding 
charge configurations near the particle-hole symmetry line. 
Filled and empty circles denote particles and holes, respec- 
tively; paired (upper and lower) solid lines represent the two 
coupled arrays, the couplings between which are illustrated 
by the dashed lines. The low-lying energy levels satisfying 
n+ (r) = are well separated by a large amount of energy (of 
the order of Eq) from those with n+(r) 7^ 0. 



(b) 

FIG. 4. Typical cotunneling processes relevant (a) near 
the particle-hole symmetry line and (b) near the maxi- 
mal-frustration line. The intermediate virtual state costs an 
energy of the order of Eq. 
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FIG. 5. Schematic phase diagram of the system. The su- 
perconducting SC phase is distinguished from the SC phase 
in that the underlying transport mechanism is the cotunneling 
process instead of the single-charge transport. 
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